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EVIDENCE OF THE BEAM PATTERN CONCEPT OF 
SUBSONIC JET NOISE EMISSION 

S. Paul Pao and Lucio Maestrello 
Langley Research Center 

SUMMARY 

This study describes the methods and results of a new approach toward positively 
identifying the apparent location and strength of the sound source in a subsonic jet. The 
Green’s function technique is applied to regenerate the three-dimensional distribution of 
the acoustic far field from measured data on a two-dimensional plane. The apparent 
sound -source location is also determined by using geometrical ray acoustics. The 
results revealed a number of features of the subsonic jet noise mechanisms. They show 
that the apparent sound source is highly localized within the jet flow and that the spon- 
taneous noise emission from such a compact source volume is in the fashion of a nar- 
row beam with random directions. These observable properties of the sound field offer 
the basis for a beam pattern concept of subsonic jet noise radiation. This study not only 
provides the immediate benefit of new insight into the mechanisms of subsonic noise radi- 
ation but also verifies a powerful and accurate technique of aerodynamic noise measure- 
ments. This technique is expected to have important applications in the development of 
jet noise reduction technology. 


INTRODUCTION 

The continual requirement for advances in jet noise reduction has provided a strong 
stimulus for research in the area of positive identification of dominant noise production 
mechanisms within the jet flow. Such work demands a substantial degree of refinement in 
both the experimental and analytical techniques. Several ways of making measurements to 
determine the sound- source location have been used: (a) direct measurement of the fluc- 

tuating velocity or density within the jet, where the sound -source strength is calculated 
according to theoretical models: (b) cross correlation between measurements of fluctuating 
aerodynamic quantities in the jet and the acoustic measurements in the far field, where 
sound- source strength distribution is determined by analyzing the cross correlation func- 
tions: and (c) acoustic measurements near the Jet, where either mechanical or analytical 
acoustic imaging techniques are used to recover the apparent sound- source distribution 
within the jet. The present study belongs to the third category. It attempts to advance the 



state of the art in this area. Some background for this study is found in references 1 
and 2. 

The acoustic field surrounding a subsonic jet is an intricate distribution of wave 
vectors. The sound-pressure fluctuation at each point near the jet is a summation of 
contributions from a large volume of the jet flow. From an experimental point of view, 
both the amplitude and phase information of the received acoustic signal must be resolved 
in order to retain sufficient information to detect the location and strength of the apparent 
sound source in the jet. Mutually independent pointwise measurements commonly chosen 
for jet noise directivity measurements cannot resolve such details. Some highly success- 
ful results in determining apparent sound-source strength distribution are reported in 
references 3 and 4 where mirror imaging techniques are employed. 

Enough data to provide a comprehensive description of the sound field can be 
obtained through the measurement of the space- time correlation function of the normal 
sound-pressure gradient over an imaginary plane near the jet. Reference 5 presents 
data obtained from such a measurement. That study computed the normal acoustic energy 
flux at the plane of measurement near the jet. Although the space- time correlation data 
contain all the necessary information to regenerate a detailed description of the sound 
field surrounding the jet, the original data alone cannot convey an intuitive and direct 
impression of the sound emission and propagation process. The basic purpose of this 
study is to retrieve meaningful numerical representations of the sound field surrounding 
this subsonic jet. 

In this study, two methods of data analysis are given. The first method is the 
direct integration approach using the Green's function technique. This approach is cap- 
able of recovering from the data quantities such as the three-dimensional distribution of 
sound-pressure intensity, the correlation functions of sound pressure in the far field, and 
the visualization of the motion of a statistical wave front. In the second method, a simple 
triangulation technique based on geometrical acoustic considerations is developed. This 
technique is relatively simple although it is capable of providing accurate information 
about the apparent sound- source location in the jet flow and the far-field sound propagation 
patterns. 

This study, by providing positive information about the subsonic jet noise emission 
pattern, is intended to contribute to research in the area of jet noise reduction. The 
methods of analysis, a description of the existing data, discussion of results, and the 
physical evidences leading to the beaming concept of subsonic jet noise radiation are given 
in the following sections. 
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SYMBOLS 


c 


D 



d 

f 


‘'O 

P 

q 

R 

Rpp 

Ryy 

^VpVp 

r 

^1 

^2 

S 

Spp 

Sqq 

®VpVp 


constant speed of sound 
diameter of jet 

differential correlation function 
distance of separation of two points 
frequency, Hz 

normal direction with respect to plane of measurement 
sound pressure 

equivalent sound- source distribution 
radial distance 

sound-pressure correlation function 
correlation of intermittency factor 
normal sound-pressure gradient correlation function 
radius measured from jet axis 
= |x - y| 

= j(x + |) - (y + v)\ 

designation of surface in three-dimensional space 

cross spectral density function of sound pressure 

cross spectral density function of apparent sound source 

cross spectral density function of normal sound-pressure gradient 




S{y),S{ri) range of integration for surface integral 
dS,dS(y),dS(r 7 ) infinitesimal surface element 
t time variable 

t' retarded time variable 

Uj jet exit velocity 

V(y),V(77) range of integration for volume integral 
W sound power 

Wj. total sound power from jet 

(X,Y,Z) local Cartesian coordinate system 

(x',y') local two-dimensional coordinate system on plane in far field 
X general three-dimensional variable in free space 

Xj distance along axis of jet 

general two- or three-dimensional variable in source region 

3 3 

dy ,d ?7 infinitesimal volume element 

Z' average direction of ray bundle 

y intermittency factor 

A finite differential quantity 

6 infinitesimal differential quantity 

T) spatial separation distance which takes both positive and negative values 

rj spatial separation vector in source region 
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6 angle measured in X,Z-plane as measured from X-axis 

^ spatial separation distance in far field 

^ spatial separation vector in far field 

differential correlation coefficients 

r,r* delay time variables 

Tj^ax delay time value at which correlation reaches maximum 

<p azimuth angle in local spherical coordinate systems 

(p* angular variable around circumference of jet 

x{/ angle in local spherical coordinate system 

00 circular frequency 


BASIC ANALYSIS 

The theory of wave propagation states that if the sound source can be described as 
a known function of space and time, the induced sound field throughout the free space can 
be computed. In some cases, however, the sound source is not known although an equiva- 
lent mathematical condition can be established. In figure 1, a distributed acoustic source 
is enclosed by a surface S in a three-dimensional homogeneous medium. If the time- 
varying normal sound-pressure gradient on this surface S induced by the source dis- 
tribution inside is known, the induced sound field throughout the free space outside of 
the surface S can be computed by using the appropriate Green’s function. Such an 
equivalent representation of the sound source is most useful in the case of subsonic jet 
noise measurement because the strength of the aerodynamic sound sources within the jet 
is not yet amenable to direct measurements. Instead of direct source measurements, 
equivalent acoustic measurements can be made on a surface within the acoustic field near 
the jet. Figure 2 shows the geometrical arrangement of such a measurement which was 
reported in reference 5. The surface S becomes simply an infinite plane separating 
the space into two halves: the inner half which contains the subsonic jet exhaust flow 
and the outer half which is source free. The acoustic data in reference 5 were obtained 
in the form of space-time correlation functions of the normal sound-pressure gradients 
on this plane. 
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The basic aim of this study is to obtain from the data meaningful numerical repre- 
sentations of the sound field surrounding this subsonic jet. Although the space-time cor- 
relation functions on the plane contain all the information to regenerate the kinematics of 
the sovmd field, the original data alone lack the ability to convey an intuitive and direct 
image of the sound emission and propagation process; consequently, further analysis is in 
order. In this section on basic analysis, the key mathematical tools for extracting useful 
information from the data are discussed. 


The Green's Function Formulation 

Derivation of the in tegral formul ae.- Analytically, the sound field in the source 
free half-space can be computed from the prescribed data of a normal sound -pressure 
gradient on the data plane if the acoustic process under consideration is deterministic. 
The Green's function of the wave equation satisfying a Neumann boundary condition on 
the plane is used in the integral formula: 


( 1 ) 


where p is the sound pressure, c is the speed of sound, x is a point in the source 
free half- space, t is the time, dS is an infinitesimal surface element on the data plane, 
y is a point on the data plane, np is the normal to the plane, and the bracket around the 
pressure- gradient function indicates that the function should be evaluated at an appro- 
priate retarded time t'. 

In this study, the sound- pressure fluctuations can be considered as a stationary 
random process. Only the statistical functional values of the sound field can be mean- 
ingful. An expression which relates the sound-pressure correlation function to the 
sound- pressure gradient correlation function can be derived: 





^vpVp(y>j?>T*) 


4^2 J Js(y),S(T7) |x - yj I X + I - (y + 2) 


dS(y) dS(rj) 


( 2 ) 



T = T - 


I'x + I - (y + 2) - X - y 
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The value Rpp(x,5,r) is the sound- pres sure correlation function as defined in the half- 
space, normal sound-pressure gradient correlation function as 

defined over the data plane, x and y are the reference points for correlation, ^ 
and 'q are the spatial separation vectors, dS(y) and dS(? 7 ) are the infinitesimal sur- 
face elements for integration on the plane corresponding to the y and the q variables, 
r is the delay time, and r* is the delay time for the source function where appropriate 
retarded time relations have been taken into account. 

The right-hand side of equation (2) is a double surface integral; each integral has an 
infinite range of integration over the entire data plane. The y -variable is fixed for the 
integration with respect to q. Therefore, the source function in this first step of integra- 
tion is simply a space- time correlation function with respect to a single reference point 
on the data plane. Such a space- time correlation function is what is commonly mea- 
sured in an experiment. A physical interpretation of this first integral can be provided. 
Equation (2) is the statistical average of the joint contribution of 3p/9nQ from an infini- 
tesimal area element at y and Op/ang from all over the data plane, while the function 
R^p^p^y,T 7 ,r*) helps to reduce the range of integration to an effectively finite domain on 
the data plane. Since all the contributions related to the sound-pressure gradient at y 
have been taken into account, this first integral represents the total acoustic contribution 
from the jet transmitted through an infinitesimal area at y. Hence, the subsequent inte- 
gration over the y -variable provides the sum of all the contributions from the acoustic 
source as transmitted through this plane. 

In this study, the available data consist of correlation functions measured at a 
limited number of reference points in the y-variable. Such data can best be used for the 
study of the local sound emission properties at various parts of the jet flow. Integration 
over the y-variable is not discussed. Since only one integral in equation (2) is evaluated, 
the equation can best be written in a differential form: 


^S(y) 


VpVp' 




. r _ 

477 ^ ^s(77)|5-y||5 + |-(y + 2 )| 


dS(77) 


It is convenient to define 


r - 


( 3 ) 


The quantity d|^p( 5,|, is called the differential correlation function and represents 
the acoustic contribution to the correlation function Rpp^x,^,T^ from a unit area in the 
neighborhood of the fixed point y on the data plane. 
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For economy of numerical computation, it is worthwhile to consider the meaning of 
this function. The function D Rpp^x,^,r) depends on seven independent variables. Con- 
sider first the case with x fixed. The restricted function then depends on four indepen- 
dent space and time variables. In general, this differential correlation function represents 
the linear dependence between sound-pressure signals at a fixed point and at a variable 
point for all values of time delay. In essence, such linear dependence can be classified 
as either time-like or space-like: 

(a) The time-like property of sound-pressure correlation between two points x 
and X + ^ states that the same acoustic pressure disturbance arrives first at one 
point X at time t and then at x + ^ at a later time t + t where r= |^|^c. The 
differential correlation function is defined as the contribution from a given unit area on 
the data plane. If x and x + ^ happen to fall along the path of propagation of the domi- 
nant acoustic energy flux coming through this unit area, a significant correlation coeffi- 
cient can be detected at very large spatial separation distances. 

(b) The space-like property refers to the correlation between the simultaneous val- 
ues of sound pressure at x and x + ^. No cause and effect relation can be established 
between these points because there is no time for any acoustic signal to travel from one 
point to another. The linear dependence of the sound pressure between x and x + ^ at 
T = 0 is a measure of the inherent similarity of the pressure field. In the present case, 
the inherent similarity results from the fact that the far-field sound pressure is induced 
by the coherent acoustic energy flux coming through a unit area of the data plane. Within 
the space and time framework, the spatial correlation pattern induced by such an apparent 
sound source propagates outwards in the free space in exactly the same manner as a wave 
front produced by a deterministic sound source. Such analytical similarity can be verified 
through the definitions of equations (1) and (2). Based on these reasons, the space-like 
correlation pattern of the differential correlation function for sound pressure in the far 
field can be defined as a statistical wave front. 

In physical situations, the general behavior of the differential correlation function repre- 
sents a mixture of the space-like and the time-like properties of the sound-pressure field. 

The equivalent sound- source distribution.- Under the Green's function formulation, 
it is possible to establish an equivalent sound-source distribution by means of a narrow- 
band spectral analysis. If an apparent sound-source distribution q(^y,t^ can be assumed 
in the jet flow, with q(y,tj being a random function, then a cross spectral density function 
can be defined for the source distribution function, which can be written as 

®qqfe-2'") = S j_„ ‘ + dr (« 
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Similarly, one can define a cross spectral density function in the far field, so that 




( 5 ) 


The functions and Sqq^y,? 7 ,a>^ can be related by means of a double integration 

formula if the simple wave equation is assumed to hold in spite of the presence of the jet 
flow. This equation can be written as 


where 



— — 





(6) 



r 


2 


X + (y + ^) 


The quantities dy^ and drj^ are infinitesimal volume elements, and both of the domains 
of integration V(y) and V{ri) are the entire three-dimensional volume over which the 
apparent source function q(y,t) is defined. The cross spectral density functions in 
equations (4) and (5) and the Green’s function in equation (6) are complex quantities. 

Equation (6) is a multidimensional linear relationship between Spp^x,^,co) and 
Sqq(y,775o;). In a practical situation, Spp^x,|,o;] rather than Sqq^y, 77 ,ct)^ is more 
often measurable. It can either be measured directly or derived indirectly from other 
measurable data. For example, Sppj'Xj^jCoj can be computed from the space-time cor- 
relation function of the normal sound-pressure gradients on the plane mentioned earlier 
in this section: 


^pp(5’i’^) - ^ JJ, 




’S(y),S(r7) r^r 2 




with 


®VpVp(y’2’‘^) ■ 27T ® ^^^^VpVp(M»'^) 

The function Sypyp^y,? 7 ,a)^ is a cross spectral density function with respect to time 
variable, and a correlation function in all of the spatial variables. Such a function can 
be measured directly in an experiment. The function Spp(x, ^,aj) evaluated at ^ = 0 
is proportional to the acoustic energy flux normal to the surface S. 
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Equation ( 6 ) can be considered as a Fredholm integral equation of the first kind if 
the roles of known and unknown functions are reversed. In the configuration shown in 
figure 2, the source region and the far field are separated by the plane S; therefore, 
r^^ and r^ are always greater than zero. It follows that the kernel function of the 
Fredholm integral equation is square integrable, and it now appears that this integral 
equation has a unique solution. This is not the case, however. The unknown function 
Sqq(y,^,w) depends on six spatial variables, with co treated as a fixed parameter. The 
left-hand side of equation ( 6 ) also contains six spatial variables. However, equation (7) 
indicates that Spp^x, 4 ,ojj is actually a function of only four independent spatial variables 
because each of the y and 77 variables on the right-hand side of equation (7) is defined 
over a two-dimensional surface. A unique solution can be constructed only if additional 
conditions are known or are imposed upon the unknown source distribution function. In 
the case of computing an equivalent distribution of apparent sources, these conditions are 
arbitrary. For example, a unique solution can be obtained from equation ( 6 ) if the source 
function is assumed to be an infinitely thin two-dimensional distribution over a conical 
surface enclosing the jet flow. The unknown quantity Sqq^y,? 7 ,co^ in equation ( 6 ) is now 
a function of four variables. Physically, this equivalent sound-source distribution is a 
two-dimensional image of the actual three-dimensional sound- source distribution within 
the jet. The sound-source strength is proportional to the acoustic energy flux density 
passing through the conical surface. Another common assumption is the consideration of 
the equivalent source as a line distribution along the axis of the jet. This assumption is a 
special case of the example just given because the axis of the jet can be considered as a 
degenerated conical surface when the apex angle reduces to zero. 

Ray Path Analysis 

Instead of analyzing the wave propagation process by means of the Green's function 
formulation, a study of the delay time relation in the space-time correlation data can pro- 
vide useful information about the geometrical structure of the acoustic field near the 
subsonic jet. A correlation function of acoustic signals between a fixed point x and any 
given point x + £ is a function of the delay time t. The geometrical configuration is 
shown in figure 3, and the definition of correlation functions can be found in equation ( 2 ). 
The delay time at which the correlation function reaches its maximum is designated at 
Tmax* time of arrival of the maximum correlation value can be interpreted as the 

difference in time for a statistical wave front emanated from the jet to touch the points 
X and X + ^. Since such a statistical wave front travels at the speed of sound in the free 
field, the direction of travel of this statistical wave front can be computed by means of a 
simple equation: 

e = cos-'(i^) (9) 
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with 


^ = |i| 

where 0 is the angle between the direction of the statistical wave front and the line 
joining the points x and x + ^. 

The derivation of equation (9) can be shown in a two-dimensional configuration with- 
out loss of generality. Let the statistical wave front be first assumed as a plane. The 
cosine relation within a right triangle leads to equation (9) and is clearly shown in fig- 
ure 3. As a second derivation, let the statistical wave front be assumed as spherical. 

The delay time relation is equivalent to a statement that the difference in distances from 
the center of the sphere to the points x and x + ^ is a constant. It is known in plane 
geometry that the locus of the center of the sphere describes a hyperbola. An asymptote to 
the hyperbola is a straight line passing through the midpoint between x and x + ^ with 
an inclination defined by the angle 9, In a general three-dimensional situation, the geo- 
metrical configuration as described in figure 3 has a rotational symmetry about an axis 
passing through the points x and x + ^. Hence, the locus of the center of the sphere is 
a hyperboloid with a conical asymptotic surface. (See sketch in fig. 3.) 

The two-dimensional configuration in figure 3 shows that the hyperbolic locus of the 
center of the sphere comes very close to the asymptote, starting a short distance from the 
points X and x + 4 . Therefore, in this report, the asymptote is termed the ray path, 
which is a practical approximation to the geometrical acoustic ray representing the local 
spherical statistical wave front. Although both the plane wave and the spherical wave 
assumptions lead to the construction of nearly identical ray paths, physical interpreta- 
tion of the results of subsonic jet noise measurement is much easier with the spherical 
assumption because it is known that the jet noise source region lies within only a finite 
distance from the plane of measurement. 

DESCRIPTION OF EXISTING DATA 
Experimental Setup 

A sketch of the experimental setup is shown in figure 4. A subsonic jet was used in 
the experiment described in reference 5 where the exit diameter was 6.25 cm, the stag- 
nation temperature was ambient, and the exit velocity was 213 m/sec. The boundary of 
the mixing flow regime was defined as the 1-percent intermittency surface around the 
jet. This boundary, based on measurements in reference 5, was a conical surface with an 
half apex angle of 11®. The acoustic measurement was made on a plane near the jet which 
made an angle of 17®20’ with the axis of the jet. The orientation and the location of the 
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plane were chosen so that the pressure fluctuations in the neighborhood of the plane were 
governed by the inverse square law, yet the entire plane remained as close to the jet as 
possible. Ten stations were chosen on the plane as the reference points for correlation 
measurements, and their locations are given in table I. The microphone placement around 
each station was arranged in a polar pattern along radial directions 45° apart. The 
typical radial separation was 4 cm, with between 75 to 120 points selected at each station. 
In connection with the present analysis, local Cartesian coordinates are defined so that the 
origin coincides with the station reference point. As stated in the preceding section, only 
the local properties of the acoustic data are analyzed, and each station is treated inde- 
pendently. The imaginary data plane serves as the X,Y-plane for the local coordinate 
system with the X,Z-plane passing through the axis of the jet. The outgoing normal to the 
plane is designated as the positive Z-axis. This local Cartesian coordinate system serves 
as either the x coordinate system in the far field or the y coordinate system on the 
data plane. Other details of the experiment can be found in reference 5. 

General Characteristics of the Data 

The pressure gradient normal to the plane was directly measured by using an 
electret pressure-gradient microphone (ref. 5). Since the linear response of this trans- 
ducer was restricted in frequency, the data are prefiltered to within a band from 250 to 
3000 Hz. The peak frequency of the jet is close to 1000 Hz; hence, approximately one and 
one-half octaves on either side of the peak of the jet noise spectrum are included in the 
data. Cross correlation functions for the prefiltered signal from pairs of microphones 
were obtained by using a real-time correlator. All of the cross correlation functions 
were taken with a 20-microsecond time delay interval with a few exceptions where a 
50-microsecond time interval was used. Typical sets of the correlation function are 
shown in figures 5, 6, and 7. 

These figures show the values of the space-time correlation function RvpVp(®»!i»'’') 
as a function of the spatial separation vector and the delay time t. The first variable in 
the correlation function is zero because the reference point of correlation is chosen as 
the origin of the local coordinate system. The result of measurement at station 1 is 
shown in figure 5. Three groups of curves are in this figure with the direction of the 
spatial separation vector t? fixed in each group. Hence, an individual curve represents 
the variation of I^vpVp(®’’^’'^) ^ function of t at a given value of the separation 

distance r] along the given direction of the group. Beyond a certain value of spatial 
separation, the value of the cross correlation function ceased to maintain any significant 
value, at which point the series of measurements would terminate. Such maximum dis- 
tance of correlation can be referred to as the range of correlation. Incidentally, the 
maximum correlation amplitude at large values of 77 occurs also at a large value of t. 
Hence, a group of correlation function curves which has a large correlation range in the 
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spatial domain also has a large correlation range in the time domain. The group of cor- 
relation curves with the spatial separation in the Y-direction in figure 5 are given only 
for positive values of rj. Measurements in reference 1 had shown that the correlation 
patterns were symmetric with respect to the X,Z -plane. For the data at station 1, the 
correlation range for each of the chosen directions is approximately the same. As mea- 
surement proceeded from one station to another, the correlation range in the Y-direction 
remained approximately the same; however, the correlation range changed significantly in 
the X-direction. Figure 6 shows three groups of correlation curves with separation dis- 
tances in the X-direction for stations 3, 5, and 7. The range of correlation increases 
progressively from the upstream stations to the downstream stations. Figure 7 shows 
the entire set of data from station 9. The great disparity of correlation range in the 
X-direction and the Y-direction can easily be noted. Significant correlation values 
remain observable for separation distances of more than 15 diameters along the positive 
X-axis. The continual elongation of the X-direction correlation range can be explained 
later in this paper when the concept of beam pattern radiation of jet noise is introduced. 

Although the above-mentioned figures contain sufficient information to describe the 
propagation of sound near the jet, these figures lack the ability to convey an intuitive visual 
impression to the reader. There is a better way of data presentation for the purpose of 
visualization. The correlation function can be given in terms of spatial contours at a 
fixed value for delay time. The data set, as obtained at station 8, is presented in the form 
of such spatial contours in figure 8. This figure contains four sets of spatial contours at 
T = -0.2, 0.0, 0.6, and 1.4 milliseconds shown from left to right in the order of increasing 
delay time. This figure shows contours which represent negative values of the correlation 
function as dotted lines. If each set of spatial contours of the correlation function is con- 
sidered as a statistical wave front (discussed in the previous section), then figure 8 
illustrates the propagation of such a statistical wave front from upstream to downstream 
along the data plane. The center of the statistical wave front touches the origin of the 
local coordinate system at r = 0. The shape of the statistical wave front appears to 
spread out and form a circular arc as the delay time approaches larger values. 

DISCUSSION OF RESULTS 
Sound-Source Location 

Construction of the ray bundle .- The method of ray path analysis is a simple way to 
find the apparent sound- source location by means of the measured space-time correlation 
function of the normal sound-pressure gradient on an imaginary plane near the subsonic 
jet. For data sets such as those shown in figures 5 and 7, one geometrical acoustic ray 
can be constructed from each time correlation function at a given spatial separation. All 
the rays at one station form a bundle. The ray bundles at stations 3 and 8 are shown in 
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figures 9 and 10, respectively. Figure 9 shows the rays reconstructed at points along the 
X-axis and the Y-axis in the local coordinates at station 3. For the bundle along the 
X-direction, the rays converge to a very small region within the jet directly opposite to 
the location of station 3. The bundle along the Y-direction also has a remarkably sharp' 
convergence to a small region. These two bundles can be fitted together to form a three- 
dimensional bundle which is illustrated by the sketch in figure 9. Only the ray bundle 
along the X-direction is shown in figure 10, In this case, the bundle converges to a some- 
what larger region in the jet. Figure 10 demonstrates the significant point that the 
location of the apparent source region within the jet is far upstream from the reference 
point at station 8. Table I shows that station 8 is located at 16,65 diameters downstream 
from the nozzle exit, while the average location of the apparent source region is approx- 
imately 8.8 diameters from the nozzle exit. 

As these figures show, there is some random deviation in the ray direction. Only 
the average direction of the ray bundle is meaningful. After constructing all the ray bun- 
dles for all the stations, parameters related to the geometrical configuration of the ray 
bundles (average direction, average location, and approximate length of the region of 
convergence, etc.) have been obtained and are summarized in table I. Figure 11 shows 
the average ray bundle direction for all the stations and the outline of the ray bundle for 
every third station. It is obvious from this figure that there is continuous variation of the 
noise emission pattern from different parts of the jet mixing flow. 

For each of the 10 stations the ray bundle converges to a relatively small region in 
the jet flow. The convergence is a strong indication that the dominant part of the acoustic 
energy reaching the plane of measurement originates from small volumes of coherent ran- 
dom sound sources, A quantitative analysis of the apparent source strength distribution 
for each of the ray bundles can be done, in principle, according to earlier discussions 
concerning the Greenes function formulation. However, such an analysis is not included 
in this study; an indirect estimate of the apparent source strength distribution is discussed 
in the following section. 

In reference 6, Maestrello and McDaid analyzed the narrow-band acoustic energy 
flux vectors near a subsonic jet according to measurements obtained in the neighborhood 
of one reference point. The reader may be interested in referring to this previous work 
for comparison with the present results. 

Comparison with the mirro^ techn ique .- The normal acoustic energy flux along the 
plane of measurement was previously computed in reference 5 based on the experimental 
values of I^vpVp(y?^^'^)* ^ further assumption of symmetry about the axis of the jet, 

the apparent source strength distribution per unit length of the jet was computed by mul- 
tiplying the acoustic energy flux density by the area of a ring segment of the conical sur- 
face which is generated by rotating the plane of measurement around the axis of the jet. 
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Figure 12 shows the results of the power per unit length normalized against the total 
power of the jet and plotted against the Strouhal number. The pealc amplitude occurs at 
fD^Uj = 0,3, and it is in good agreement with the commonly measured peak value of the 
far-field jet noise spectrum. In the spatial coordinate, the variation of the apparent 
source strength is very gradual except for the initial few diameters. The peak apparent 
source strength per unit length of the jet occurs at approximately 9 diameters downstream. 
However, such apparent source strength distribution is obtained according to the normal 
acoustic energy flux at the plane of acoustic measurement. If the vectorial property of 
the acoustic energy flux is considered, the location of the apparent source should be shifted 
from its location along the plane to its apparent location along the axis of the jet. The 
shift can best be made by following the average direction of the ray bundle. The shifted 
apparent source strength distribution has an important advantage: it is independent of the 
position of the plane of acoustic measurement. For example, if a shift is applied to the 
source distribution of figure 12, the point of maximum strength is located at approximately 
6.6 diameters downstream from the nozzle exit of the jet. 

It is interesting to compare the result of this study with the mirror technique as 
recently described by Grosche, Jones, and Wilhold (ref. 3) and by Laufer, Schlinker, and 
Kaplan (refo 4). The mirror is a high-gain, direction- sensitive device which has a large 
aperture to receive the acoustic signal from the jet. The measured signal represents the 
radiated sound from a small volume of the flow which is located at one focal point of the 
mirror. On the other hand, the technique described here for measuring the space-time 
correlation function of the normal sound-pressure gradient on a plane can record the 
phase and amplitude of the acoustic signal as emitted from the jet over a large area on 
the plane of measurement. The analysis of such data by means of either the delay time 
technique or the Green*s function formulation can reconstruct the acoustic image of the 
apparent sound source. Hence, the results produced by the present technique and those 
produced by the mirror technique are similar in nature. The relative advantages and 
drawbacks for these two techniques should be an important subject for further studies. 

Add itional evidences of sound- source localization .- Some supplemental experimen- 
tal data are available to support the hypothesis that a coherent volume of random sound 
sources within the jet is highly localized. These data include measurements made inside 
and immediately outside the shear layer of the jet. Within the shear layer, a set of hot- 
wire anemometers with standard signal conditioning devices to discriminate between 
turbulent and nonturbulent regions was used; outside the shear layer a pair of microphones 
were chosen for pressure measurements. The hot wires were placed around the circum- 
ference of the jet at the radial distance where the intermittency value was y = 0.5. The 
radial component of the velocity field thus obtained can be considered as the radial com- 
ponent of the local vorticity according to Corrsin and Kistler (ref. 7). Hence, the resul- 
tant spatial correlation R^0*,Xjy^D^ can be interpreted as the spatial correlation of the 
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radial vorticity component about a radius where the mean intermittency signal is 50 per- 
cent turbulent. Such data were taken at two stations downstream of the nozzle at 
Xjy^D = 4 and Xj^D = 8; the results are shown in figure 13. 

The correlation is of appreciable amplitude only over azimuth angles of about 10° 
in either the +0* or -(p* directions. Although this result does not describe a com- 
plete source element in the jet, it does support the view that the true sources are local- 
ized and uncorrelated over large azimuthal angular separations. This result again 
supports the view that localized sources in the flow produce the sound field. 

Immediately outside the jet, the pressure fluctuations can be thought of as associated 
with the vorticity in the flow. Microphone correlation measurements over cp* for 
Xjy/D = 2 and Xjy'D = 4 are shown in figure 14. At 90°, the correlation level drops to 
0.4 of the peak v^ue, and becomes very small beyond 120°. It is not surprising that the 
pressure correlates over a greater angle because the microphones may receive signals 
from sources at a distance, where the propagation effects naturally broaden the correla- 
tion distance. One interesting aspect of the pressure measurement is that the peak of 
correlation for all angles occurs at zero time delay. This evidence tends to rule out the 
possibility of dominant spiral modes in the flow. In a recent paper by Michalke and Fuche 
(ref. 8), similar observations and conclusions were reported. Additional evidence that the 
pressure about the jet is not fully correlated around the circumferential direction was 
reported by Liu, Maestrello, and Gunzburger (ref. 9). 


Statistical Nature of the Sound Field 

The fa r-field distribution of the me^ square sound pressure.- A highly descriptive 
numerical result is the distribution of D Rpp(x,0,0)j in the far field. The result is the 
mean square pressure distribution as induced by the acoustic energy flux coming through 
a unit area on the data plane at a reference station. The contour plots of D|Rpp(|x,0,0jj 
in the X,Z -plane for all the odd number stations are shown in figure 15. The values of the 
contour are given in physical units, and the coordinates are given in terms of X/D and 
Z/D. The general direction of wave propagation agrees very well with the ray path 
analysis. If the sound wave passing through each station is assumed to originate from the 
tip of the corresponding ray bundle, the distribution of D Rpp(x,0,0) follows the inverse 
square law rather accurately. 

In a previous paper describing the initial results of this study (ref. 1), some negative 
values of DjRpp^x,0,0^ were noted in the computed results following directly from the 
application of equation (3). An example is shown in figure 16. There are two possible 
causes for the presence of such negative values; 


16 



(a) The value of D|^Rpp^x,0,0)J is a summation over all the data at each station. 
Theoretically, the integr^ should always have a positive value. If the data are imperfect, 
the numerical integrations need not be positive. 

(b) The acoustic signal is taken within a frequency band of approximately one decade. 
The phase effect of such data may have resulted in the appearance of the negative loops. 


Further investigation into the numerical details revealed that the actual value of 
Q shifted slightly in time. The average time shift is approximately 
120 microseconds, while the actual value depends on the radial distance of the point under 
consideration. A typical example is shown in figure 17. Such a time shift is very prob- 
ably a direct result of data imperfection compounded by the strong modulation of the cor- 
relation function which was caused by the band limitation of the original acoustic signal. 


The direction of flow of the dominant portion of acoustic energy can be determined 
roughly from the distribution of DjRpp^x,0,0) . On a point-by-point basis, the local wave 
propagation direction can be determined rigorously by computing the maximum of the 
local pressure-gradient correlation function. An integral formula for such computation 
can be readily derived from equations (1) and (2) by differentiating the entire equation 
with respect to the direction of the desired component of pressure gradient. However, the 
formula would involve a second time derivative of the function RvpVp(y»^>'^) which is a 
measured quantity. Hence, the numerical accuracy of such computation is questionable. 


The correl ation relations of sound pressure in the far field .- The cross correlation 
function of sound pressure provides a quantitative measure of the linear dependence of the 
acoustic signal at different points. The normalized form of the differential cross cor- 
relation function provides a convenient way to study the statistical structure of the induced 
sound-pressure field. The normalization can take two different forms: 



-1 


(10) 


and 


(D Rpp|^x,0,0jj DjPtpp^^x + ^ D Rpp(X,^,Tj 


- 1/2 


( 11 ) 


The first definition is a normalization of the differential correlation amplitude with respect 
to its maximum value at the reference point. The second definition has taken into account 
the local maximum amplitude at both points as involved in the cross correlation. The 
second definition is preferred from a physical standpoint. 


A typical set of reconstructed differential correlation functions are shown in fig- 
ure 18, which are based on the data at station 5. The reference point and the variable 
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point are confined to the X,Z -plane. The moving spatial point is taken along a line perpen- 
dicular to the direction as defined by the origin of the local coordinate and the reference 
point. The geometrical configuration is shown in a sketch in figure 19. 

Since the wave propagation process becomes self- similar in the far field along the 
radial 'direction in a spherical coordinate system, the spatial distribution of the differential 
correlation coefficient can best be described in such a coordinate system (fig. 19). Con- 
sider the case of a reconstructed differential correlation function where both the reference 
point and the moving point are confined to a spherical coordinate surface. The spatial 
separation can then be described in the angular coordinates ^ and cp. In many of the 
computations, the correlation coefficient is computed along lines tangent to the sphere 
at the reference point. The correspondence between the angular coordinate cp and the 
linear separation distance can be given simply by: 

tan (p = ^ d = I ^ I (12) 

An inverse square law correction factor must be applied to the differential correlation 
coefficient Pq^x,^,t) when the description is transferred from the tangent to the spherical 
surface: 




(13) 


The value of remains unchanged for the same coordinate transfer. 

Values of p {x,cp,T) and p_^{x,(p,T^ are computed at all of the stations. The 
reference point of correlation is typically chosen at a distance 100 diameters away from 
the origin and along the average outgoing direction of the ray bundle at this station. For 
each given value of cp the correlation coefficient reaches a maximum at a certain delay 
time T . The envelope of such maximum values provides a measure of the width in 
the sound field over which the sound-pressure field fluctuations remain linearly dependent. 
Note that the boundary of correlation is now expressed in angular coordinates. As long as 
the reference point is chosen along the average ray bundle direction and the radius is very 
large, the values of PQ[^,4>Jrna^ P max) invariant. Hence, the domain 

of correlated sound propagation is in the pattern of a narrow beam as confined within a 
conical surface. A detailed discussion of the beam pattern concept is given in the follow- 
ing section. 


Both PQ(x,</),Tj^ax) Pl(5>'?^j'^max) rapidly decreasing functions of (p. 

The actual dependence varies from station to station. The results are summarized in 
figures 20 and 21. In these figures, the average is shown in circles and the maximum 
range of variation of the differential correlation coefficients is shown by means of a line 
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interval. A curve for cos"^ <p is shown in both figures for comparison. This particular 
function was chosen because it represents the directivity of an idealized quadrupole. 

Figure 20 shows that the average shape of pQ(x,0,Tjpg^) is substantially narrower than 
cos^ <p. The angular span between the points where the differential correlation coefficient 
reaches 0.5 is approximately 25°. Figure 21 shows that the value of Pj(x,<|>,Tjyjax) 
initially follows the cos'^ curve, but then it goes above this curve for larger values 
of (p. The angular span between the points = 0.5 is approximately 35°. 

A particular set of computations based on the data at station 8 was obtained in order 
to show the statistical structure of the pressure field on a plane which is perpendicular 
to the average ray bundle direction. At station 8, the average direction of the ray bundle 
is 33° from the data plane. The reference point of correlation is chosen as 100 diameters 
away from the origin. The plane in which 4 varies is tangent to the coordinate sphere 
passing through the reference point. Local coordinates are given on this plane in order 
to fix its orientation with respect to the Cartesian coordinate system at station 8. The 
X'-axis lies within the X,Z-plane, and it points in the general downstream direction of 
the jet. The Y’-axis is normal to the X,Z-plane. In figure 22, two groups of curves are 
shown. At each point in this plane, the maximum correlation amplitude arrives at a 
specific delay time Tj^ax' contours of constant time of arrival are given as the 

dashed lines in figure 22. Note that this group of contours forms a set of almost perfect 
concentric circles around the origin. This formation is typical of the traces produced by 
the passage of a spherical wave front through the plane. In particular, the direction of 
propagation is perpendicular to this x',Y*-plane. The other set of contours, as shown in 
solid lines, denotes the distribution of the maximum value of the differential correlation 
function d| Rpp(x,^,T)] . This second set of contours shows that the amplitude of corre- 
lation is directionally dependent. In general, the correlation amplitude decreases 
rapidly with increasing spatial separation in the x'-direction; such a decreasing trend 
is somewhat slower in other directions on the plane. 

Figure 23 shows another view of D Rpp(x,§,r) in the far field. The differential 
correlation function is computed with a fixed reference point A and a fixed delay time 
T= 1.0 millisecond. The value of 4 is restricted to the X,Z -plane within an annular sec- 
tor as shown by the heavy dashed line in figure 23. As previously discussed, the spatial 
correlation contours can be considered as a statistical wave front. In this case, the sta- 
tistical wave front is induced by the acoustic energy flux from a unit area on the data plane. 
Since the point A is chosen as the reference point, this statistical wave front reaches A at 
T = 0. It moves away into the far field as t increases. The arc F-F would indicate 
the position of the correlation contour pattern if it is computed at a delay time somewhat 
greater than 1.0 millisecond. The average position of this correlation contour pattern 
is approximately 35 diameters away from station 8. Remarkably, its shape is nearly 
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circular within a short distance from the jet flow. This is a strong indication that coher- 
ent sound sources in the jet are highly localized. 

The beam pattern concept .- The above analysis allows a number of new features of 
jet noise mechanisms to be confirmed. First of all, the location of the apparent sound 
sources can be determined accurately. The source region within the jet responsible for 
the emission of correlated sound in the far field appears to be highly localized. Further- 
more, the sound pressure in the far field is correlated within a narrow conical domain. 
This domain can be referred to as a beam pattern. This use of the term "beam pattern" 
differs from conventional usage where the width of the beam is normally defined accord- 
ing to the acoustic intensity. In the present case, the statistical beam width is defined 
according to the value of the differential correlation coefficient across the beam. 

Both the correlation measurements and the statistical properties derived from the 
data provide strong evidence that some order does exist in the random acoustic field near 
a subsonic jet. Since assumptions drastically contrary to existing theories of jet noise 
are unnecessary, the following physical model of subsonic jet noise radiation can be 
postulated. 

Jet noise is radiated from localized volumes of correlated turbulence within the jet 
flow. Within a given short period of time, the coherent noise radiation is pointed mainly 
in one direction in the fashion of a narrow beam. However, the direction of the beam 
from such a radiating volume is random in time. Therefore, the commonly measured 
smooth directivity pattern is an envelope of a large number of such narrow beams. 

The time randomness of the beam pattern is simply a reminder of known properties 
of turbulence. Since the correlation functions are measured by means of a time averaging 
process, the time randomness property is not distinguishable in experimental measure- 
ments. It is important to point out such a possibility since the common emphasis in the 
literature for the purpose of explaining the structure of the acoustic far field is placed on 
the spatial randomness of the distribution of turbulence. With the assumption of time 
randomness, the size of the spatially coherent sound- source volume can either be large 
or small while the overall directivity as produced by such sound sources can remain 
compatible with known experimental evidences. 

Such a description of jet noise radiation mechanism is also consistent with the con- 
cept of a turbulent quadrupole. In jet noise theories, the overall contribution from all 
components of the quadrupole tensor to a point in the far field is equivalent to a single 
longitudinal quadrupole pointing in the direction of the receiver. As the data analysis 
has shown, the directivity of the narrow beam pattern is indeed close to the shape of 
cos^ 0. 

An immediate application of the beam pattern concept is the provision of a physical 
interpretation of the properties of the original space-time correlation data of the normal 
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sound-pressure gradient on a plane near the jet. For the upstream stations, the direction 
of the beam pattern is almost perpendicular to the plane. Hence, the correlation patterns 
on the plane represent the relatively short space-like correlation range across the narrow 
beam. For the downstream stations, the direction of the beam makes a very small angle 
with the surface of the plane. The statistical wave front within the narrow beam propa- 
gates from one end of the plane to another sweeping out into a much elongated correlation 
pattern along the X-direction. 

The beam pattern concept together with the methods of analysis described in this 
study should aid the development of techniques to observe further details of the sound- 
source structure within the jet. 


CONCLUDING REMARKS 

By using a new approach in analysis, the information retrieved from measured 
space-time correlation functions of the normal sound-pressure gradient on an imaginary 
plane near a subsonic jet has suggested the following conclusions: 

1. The apparent origin of coherent sound emission within the jet flow is localized. 

2 . The apparent sound-source location can be determined accurately by reconstruct- 
ing a geometrical acoustic ray bundle at each station of acoustic measurements. Each of 
the ray bundles converges to a specific region in the jet flow. The average direction of 
the ray path varies continuously from station to station. The most important conclusions 
here are: (a) in the near field, the acoustic energy flux is not normal to the jet boundary, 
and (b) the apparent sound-source locations in the jet may lie at a point significantly far 
upstream from the reference point where the acoustic measurements were made. 

3. Within the reconstructed sound-pressure field far away from the plane of measure- 
ment, the sound pressure is correlated only between points confined within a conical region 
surrounding the average direction along which the dominant portion of the acoustic energy 
flux is transmitted. The reconstructed sound-pressure distribution mentioned above is 
induced by the acoustic energy flux passing through a unit area on the data plane. The 
distribution of the correlation coefficient within this conical region is very similar to 

one produced by a longitudinal quadrupole pointing along the axis of the conical region, 

4. The definition of the reconstructed sound-pressure field permits the spatial con- 
tours of the reconstructed differential correlation function to be considered as a statistical 
wave front. The computed motion of such statistical wave fronts can provide valuable 
insight into the process of sound propagation near the jet. 

5. A physical model of subsonic jet noise radiation can be postulated according to 
the observable mechanisms of jet noise emission as obtained in this study. Jet noise is 
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radiated from localized volumes of correlated turbulence within the jet flow. Within a 
given short period of time, the coherent noise radiation is pointed mainly in one direction 
in the fashion of a narrow beam. However, the direction of the beam from such a radiating 
volume is random in time. Therefore, the commonly measured smooth directivity pat- 
tern i^ an envelope of a large number of such narrow beams. 

Since the data of the space-time correlation measurement cover 10 stations located 
between 0 to 20 diameters downstream from the nozzle exit, continuous variation of the 
local noise emission pattern along the jet is clearly visible from the reconstructed infor- 
mation concerning the sound field. The present technique of analysis has proven to be a 
dependable practical method for comprehensive studies of sound emission mechanisms in 
the noise production region of the jet. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
December 5, 1975 
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TABLE I.- STATION COORDINATES AND THE LOCATION 
OF APPARENT SOURCES 
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Figure 4.- Setup of pressure-gradient microphone measurements near subsonic jet 
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Figure 5.- Space-time correlation functions of the sound-pressure gradient at station 1. 
Spatial separation 77 is measured along (a) Y-axis; (b) 45°-direction; and (c) X-axis. 
Units are measured in 10"^^ N^/cm^, 
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Figure 6.- Evolution of space-time correlation function from station 3 to station 7 
The spatial separation is measured along the X-axis. Units are measured in 
IQ-IO N^/cm®. 
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Figure 7,- Space-time correlation functions of sound-pressure gradient at station 9. Spatial separation 
measured along (a) Y-axis; (b) 45° ray; and (c) X-axis. Units are measured in 10"^*^ N^/cm®. 









Figure 10.- Ray bundle reconstruction at station 8 as shown on meridian plane of jet. 










Figure 13.- Measurement of intermittency correlation coefficient around jet. 
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Figure 14.- Measurement of pressure correlation 
coefficient around jet. 
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(d) Station 7. (e) Station 9. 

Figure 15.- Reconstructed far -field sound-pressure intensity distribution transmission though a unit 
area at stations 1, 3, 5, 7, and 9. Time shift correction was included. The unit on the contours is 






X 10"^ 



Figure 18.- Typical set of reconstructed sound-pressure correlation function 
DfRpp(x, Reference point x was chosen to be a point lOOD away 

from station” 5 along a ray 67° from data plane. Units are measured in 
10"1® N^/cm^. 
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Figure 20.- Values of across beam pattern. 





Figure 22.- Contours of amplitude and time of arrival of maximum sound -pressure 
correlation on plane normal to ray bundle axis. This computation is based on data 
from station 8. Units are measured in 10“^® N^/cm®. 
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